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Abstract 

Derealization problem for a two-dimensional non-interacting electron system 
is studied under a random magnetic field. With the presence of a random 
magnetic field, the Hall conductance carried by each eigenstate can become 
nonzero and quantized in units of e 2 /h. Extended states are characterized by 
nonzero Hall conductance, and by studying finite-size scaling of the density of 
extended states, an insulator-metal phase transition is revealed. The metallic 
phase is found at the center of energy band which is separated from the 
localized states at the band tails by critical energies ±E C . Both localization 
exponent and the critical energy E c are shown to be dependent on the strength 
of random magnetic field. 
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The Anderson localization theory 1 ' 2 predicts that all states in a two-dimensional (2D) 
electron system are localized in the absence of a magnetic field. The quantum Hall effect 
(QHE) system is a first example of 2D systems which show the existence of truly extended 
states. 3 ' 4 In this latter case, the presence of a magnetic field breaks the time-reversal symme- 
try and destroys constructive interference of the backward scattering 2 so that it is possible 
for electrons to propagate forwardly. 

Recently, an intensive attention has been attracted to the derealization problem in a 
2D random-magnetic-field system. This problem is closely related to the half-filled QHE 
system 5 ' 6 as well as the gauge-field description 7 ' 8 of the high-T c superconductivity prob- 
lem. However, despite a lot of numerical and theoretical efforts, the issue of derealization 
still remains controversial. Theoretically, Zhang and Arovas 9 have recently argued that the 
field-theory description, which corresponds to a non-linear sigma-model of the unitary class 
without a topological term due to zero average of magnetic field, should have a term rep- 
resenting a long-range logarithmic interaction of the topological density ( due to the local 
magnetic field). This singular term may lead to a phase transition from localized state to 
extended one. But it is contradictory to the conclusion that all the states are localized ob- 
tained by Aronov, Mirlin and Wolfle 10 in a similar approach. Earlier numerical works 11-13 
also have given conflicting results. Recently, with a larger sample size, Liu et al. 1A have 
found a scaling behavior of the localization length near the energy band tail, which can be 
extrapolated to give a insulator-metal transition energy E c . Nevertheless, a metallic phase 
has not been directly confirmed since no scaling behavior has been found there. In the 
possible metallic region, an even larger sample size may be needed in order to distinguish 
whether the states are really extended or very weakly localized 11 with the localization length 
much longer than the sample size. 

Thus it would be desirable to study this derealization problem from an alternative nu- 
merical method which directly probes topological properties of system with less finite-size 
effect. Thouless and co-workers 15 and others 16 ' 17 have found that derealization property 
of a wavefunction in the presence of magnetic field can be well characterized by its associ- 
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ated quantized Hall conductance. Nodes of an eigenstate wavefunction with nonzero Hall 
conductance can move freely and cover the whole real space when one continuously change 
the boundary condition. 17 Such a covering of real space by the nodes has been related to 
a topological invariant integer (known as First Chern number ), which is identical to the 
quantized number of the Hall conductance (in unit of e 2 /h). Thus a nonzero Hall conduc- 
tance describes the extensiveness of a wavefunction. In contrast, a zero-Hall conductance 
state will always be localized in 2D with the presence of weak-impurities (Anderson localiza- 
tion). In the QHE system, Huo and Bhatt 18 have calculated the boundary-phase-averaged 
Hall conductance for each eigenstate of a noninteracting electron system in the presence of 
strong magnetic field, and extrapolated the density of extended states (with nonzero Hall 
conductance) to the thermodynamic limit ( sample size varying from 8 to 128). They have 
found that all extended states collapse to a single energy E c at the center of the Landau 
band, with a localization length £ oc 1/\E — E c \ u and localization exponent v = 2.4 in 
agreement with previous known results. 

In the present random magnetic field case, the total Hall conductance on average has to be 
zero. But one still finds nonzero quantized Hall-conductances for eigenstates at each random 
flux configuration. Due to the general relation between a nonzero Hall conductance and 
derealization of the corresponding wavefunction, 15 one can use this topological quantity to 
characterize delocalized states. Similar point of view also lies in the heart of the field-theory 
approach of Zhang and Arovas. 9 In this Letter, we shall use this topological property in our 
numerical approach. By studying the sample-size dependence of the density of extended 
states which are states with nonzero Hall conductance, an insulator-metal phase transition 
will be revealed. The extended states are found near the center of energy band, and the 
states at the band tail are all localized with both localization exponent and transition energy 
±E C depending on the random magnetic field strength. 

We consider a tight-binding lattice model of noninteracting electrons under a random 
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magnetic field. The Hamiltonian is defined as follows: 

H = - £ e^c+c,+E^c+Q (1) 

<ij> i 

Here cf is a fermionic creation operator, with < ij > referring to two nearest neighboring 
sites. A magnetic flux per plaquette is given as <fi = J2n a ij, where the summation runs over 
four links around a plaquette. We study the case in which for each plaqutte is randomly 
distributed between — h^n and hoir, and Wi is also a random potential with strength < w. 
For simplicity, we assume no correlations among different plaquettes for (f) and different sites 
for Wi (white noise limit). The total flux for each random configuration is always chosen 
to be zero. The finite system is diagonalized under the generalized boundary condition 
\^(i + Lj) >= e l6j \^(i) > (j = l,2 represent x and y direction respectively) with lattice width 
Li = L 2 = L, and a total number of lattice sites (sample size) is Af = L x L (the lattice 
constant is chosen to be the unit). 

The Hall conductance can be calculated by using the Kubo formula. One may relate a 
Hall conductance to each eigenstate \m>: 

(m) _ ie 2 h ^ < m\p x \n >< n\p y \m > - < m\p y \n >< n\p x \m > 

where p is the velocity operator defined as p T = iJ2i( c t+ T c i emi+T ' i ~ cf Ci +T e^ mi+T ' i ) with 
r = x or y. The total Hall conductance for the system is given by a H = J2 Sm < £F a ^ zero 
temperature, with Bp as the Fermi energy. a H will always be zero on average in the case of 
a random magnetic field. However, can be nonzero for each random-flux configuration 
because of the breaking of time-reversal symmetry. As pointed out before, a state with 
nonzero (quantized) Hall conductance represents an extended state in the thermodynamic 
limit, whereas a zero-Hall conductance state should always be localized in 2D. 

Direct calculation of matrix elements in formula (2) is time consuming. We can make 
a unitary transformation |$ >= e ~ ieiX ^ L e' t92y ^ L \^ > such that the new state |$ > satisfies 
a periodical boundary condition. The Hamiltonian (1) is transformed in terms of au +T — > 
au +T + 9 T /L with r = l(x) or 2 (y). Then the boundary-phase averaged conductivity can 
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be related to partial derivatives of the wave function in the following form: 15 

^|>e(w^-^mm) <3) 

where the closed-path of the integral is along the boundary of a unit cell < 9i,9 2 < 2n. 
affl in (3) can be shown 15 ' 16 to be quantized in unit of e 2 /h. Here $ is required to be an 
analytic wave function in 2D #-space. Starting from the wavefunction $(0, 0) defined at 
one corner of the boundary in #-space, the phase of the wavefunction $ can be uniquely 
determined 15 by a process of parallel translation, first along 6\ axis and then along 6> 2 axis 
as: 

p. ((M )^.0, ,4a) 

X>*(eiA)^|^ = o. (46) 

Numerically we have diagonalized the Hamiltonian with boundary angle varying in whole 
2n x 2n phase space for each given random flux and potential configuration. At each step, 9j 
may only change by a very small value such that d<&/dQ can be well approximated by A$/A# 
(Usually A9 < 27r/100 which is adjustable in our numerical calculation to give a reliable 
result). By constructing a wavefunction satisfying conditions (4), the Hall conductance 
averaged over boundary angle is determined in terms of (3) for each eigenstate. An eigenstate 
with nonzero Hall conductance is defined as an extended state, and the corresponding density 
of states p ex t{s, A/") is obtained as a function of energy e and sample size J\f which is averaged 
over random flux-potential configurations (200 — 2000 random configurations depending on 
sample size). 

The total density of states p(e,Af) and the extended one p ext (s,N) are obtained as a 
function of energy e and lattice size Af ( A/=16, 36,64,100 and 144). The total density of 
states does not change much with lattice size, but the extended part of the density of states 
shows distinctive behaviors at different energy regions separated by critical energies ±E C . 
The ratio p ext (e,Af)/p(e,Af) is presented in Fig. 1 around — E c (with random magnetic field 
and impurity strengths chosen as ho = 0.6n and w = 1.0, respectively). All the curves in 
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Fig. 1 cross at a fixed-point e = —E c , which is independent of the lattice size. At energy 
e < —E c , the extend-state density is continuously suppressed, and can be extrapolated 
down to zero as lattice size becomes infinity (see below). On the other hand, in the regime 
— E c < e < E c , pext/p monotonically increases with lattice size and eventually saturates. 
Therefore, Fig. 1 clearly shows a metal-insulator transition at critical energies ±E C (the 
curves in Fig. 1 are symmetric about e — 0). 

Let us consider in detail the localization at the band tail e < —E c . One may define 
two quantities characterizing the localization effect: a ratio Rq of the number of extended 
states divided by the total number of states at energy e < —E c region, and a mean width 
AE of the extended states in such a regime, both of which presumably will approach to 
zero in the thermodynamic limit. Here R = JZ^ C p ex t(s)dej JZ^ C p(e)de and (AE) 2 = 
Jl^ c [£ — (— E c )] 2 p ext (e)de/ JZ^ C p ex t(^)de. R and AE versus the sample size are shown in 
Fig. 2 in a log-log plot at ho = 0.6 and w = 1.0. The data follow two parallel straight 
lines nicely, suggesting the following power-law behavior: R ~ N~ x and AE ~ Af~ x , with 
x = 0.2 ± 0.02. Such a scaling law ensures the absence of the extended states outside the 
energy range (—E c , E c ) in the thermodynamic limit. In the localized region, localization 
length is a characteristic length scale ( scaling parameter 19 ), and for a finite-size sample with 
a width L the states with a localization length £ > L should appear as extended ones. If 
the localization length goes as l/\e — (—E c )\ v when e — > —E c , one expects \/(AE) v ~ L, 
or AE ~ M- 1/2v (Ref. 18). One also has R oc AEp(-E c ,Af) oc AE. So the finite-size 
scalings of Rq and AE in Fig. 2 imply a power-law behavior of the localization length £: 
£ oc l/\e - (-E c )\ u with v = l/2x = 2.5 ± 0.3. 

At — E c < e < E c , a monotonic increase of p ex t/p with sample size is manifestly metal- 
lic behavior. It is consistent with the behavior of d(^ L / L) / dL > ( is the so-called 
decay length and (,l/L describes the extensiveness of the system) found in the metal- 
lic region of 3D system and 2D system with spin-orbit interaction (symplectic class). 19,20 
In the present approach, the quantity p ext directly characterizes extended states and can 
be extrapolated to a finite value at large sample size limit. One may also define a ratio 
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N ext /N = Pextde/ J_e c pde, namely the total number of extended states divided by the 
total number of states within {—E c , E c ). N ext /N is found to saturate to a finite value 
R c ~ 0.68 in the following manner: N ext /N — R c oc —N~ y (y ~ 0.3), which is shown in Fig. 
3 by a log-log plot (ho = 0.6, w = 1.0). A finite R c in the thermodynamic limit is a direct 
evidence for derealization. 

Very similar behaviors have also been obtained at other random-flux strengths: h = 0.4 
and 0.5 (with w = 1.0). Correspondingly, E c = ±3.7 and ±3.5, while v = 1.25 ± 0.3 and 
1.75 ± 0.3, respectively. The results suggest a non-universal localization exponent u, which 
increases with h and is consistent with a larger v (~ 4.5) obtained at ho > 0.7 in Ref. 
14. The reduction of the metallic region (—E c , E c ) indicates that the extended states are 
less favorable at larger h . With the increase of h we find that p ext becomes less sensitive 
to the sample size. When h > 0.7, dp ext /dN is relatively small around the center of the 
energy band and a larger lattice size is needed in order to get conclusive results about 
derealization. Since we study the density of states for the extended states characterized by 
the topological properties of wavefunctions (Hall conductances) , finite-size effect is expected 
to be less important here in comparison with other approaches. The existence of the fixed- 
points ±E C , which are independent of lattice size, as well as the finite-size scalings on two 
sides of ±E C indeed support this expectation. 

In conclusion, we have unambiguously demonstrated the existence of a derealization 
region for a non-interacting 2D electron system under a random magnetic field. Critical 
energy E c of metal-insulator transition has been determined. Two branches of finite-size 
scaling are found in both metallic and localized regions, and the results are extrapolated 
to the thermodynamic limit. The localization length at the band tail (e < —E c , e > E c ) 
behaves like £ ~ l/\e ± E c \ v ', with both E c and v varying with the strength of random 
magnetic field. 
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FIGURE CAPTION 



Fig. 1. The ratio of the density of extended states over the total density of states is 
plotted as a function of energy (around — E c = —3.3) for different sample size Af. 

Fig. 2. R , the number of extended states divided by the total number of states at ( 
e < —Ec), and AE, the mean width of the band of extended states in the same energy 
region, vs. sample size Af on a log-log scale. E is a mean width of the band of total states 
at ( e < —E c ). All the data are fit into two parallel straight lines. 

Fig. 3. N ext /N Q , the number of extended states over the total number of states within 
energy region (-E c , E c ), vs. sample size Af. The ratio is extrapolated to a finite number 
R c = 0.68 in the thermodynamic limit, and solid line is a best fit to the data. 
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